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We study the t-t'-t"-J model with parameters chosen to model an electron-doped high temperature 
superconductor. The model with one, two and four charge carriers is solved on a 32-site lattice using 
exact diagonalization. Our results demonstrate that at doping levels up to 1 = 0.125 the model 
possesses robust antiferromagnetic correlation. When doped with one charge carrier, the ground 
state has momenta (±7r, 0) and (0, ±7r). On further doping, charge carriers are unbound and the 
momentum distribution function can be constructed from that of the single-carrier ground state. 
The Fermi surface resembles that of small pockets at single charge carrier ground state momenta, 
which is the expected result in a lightly doped antiferromagnet. This feature persists upon doping 
up to the largest doping level we achieved. We therefore do not observe the Fermi surface changing 
shape at doping levels up to 0.125. 



I. INTRODUCTION 

It is well known that electron-doped high materials 
have very different properties compared to hole-doped 
ones. Like hole-doped materials, their undoped parent 
compounds are insulators with antiferromagnetic spin or- 
der. But the electron-doped cuprate Nd2-a;Cej;Cu04 
remains an antiferromagnetic insulator up to doping 
level X = 0.13 whereas in the hole-doped cuprate 
La2_a:Sra;Cu04 a relative small doping level of a; ^ 0.02 
is enough to destroy its antiferromagnetic correlationii 
Theoretically it has been postulated that this asymme- 
try in properties of electron- and hole-doped materials 
can be modeled by adding intra-sublattice hopping terms 
to the t-J model. ^'^ Compared to the nearest neigh- 
bor hopping motion in the t-J model, intra-sublattice 
hoppings do not frustrate the spin background. Con- 
sequently the antiferromagnetic order of the undoped 
system is better preserved upon doping. Within this 
context, various theoretical and numerical studies have 
confirmed that the electron-doped model has robust an- 
tiferromagnetic order. In addition, appropriate intra- 
sublattice hopping terms shift single-carrier ground state 
momenta from (±7r/2, ±7r/2) in the t-J model to (tt, 0) 
and its equivalent points. This means that in a lightly 
doped system small charge carrier pockets will form at 
(tt, 0) and its equivalent points in the first Brillouin zone. 
This agrees with angle-resolved photoemission (ARPES) 
experimenliii^ on Nd2-a;Cea;Cu04. Various properties of 
the electron-doped modelSi^iSiiS including its electronic 
states, spin dynamics, and Fermi surface evolution have 
been worked out with emphasis on comparison with ex- 
perimental results .JA 

In this paper we are interested in the interaction among 
charge carriers doped into the parent compound. For this 
reason we conduct a systematic study on the electron- 
doped model with a few charge carriers using exact di- 
agonalization. In this approach, larger lattices are al- 
ways preferred in order to minimize finite-size effects. 
Furthermore, in order to study Fermi surface evolution 
the lattice must have allowed k points along the antifer- 
romagnetic Brillouin zone (AFBZ) boundary, i.e., from 



(7r,0) to (0,7r). Square lattices having only 20 or 26 
sites with periodic boundary conditions do not have this 
property. The 32-site lattice is the next available square 
lattice that has this property and on which calculations 
are still manageable. The t-J model with up to four 
charge carriers on this lattice has been studied in detail 
using exact diagonalizationiiSiiiiiii^ But previous calcu- 
lations on the electron-doped model on this lattice have 
been limited to two charge carriers only.^^ Since antifer- 
romagnetic order in electron-doped cuprates is robust, 
we need more charge carriers to make the antiferromag- 
netic phase unstable. In this paper we report results for 
the electron-doped model with one, two, and four charge 
carriers on a 32-site lattice, covering doping levels up to 
X = 0.125. 

Our paper is organized as follows. We first define the 
model in section ^1 In section Hill we look at quasipar- 
ticle properties of a single charge carrier doped into the 
system. Next we consider systems with multiple charge 
carriers and study the binding energies in section Hvl and 
the real space charge carrier correlation in section 
They provide the first evidence that charge carriers are 
unbound at these doping levels. Section IVII deals with 
the momentum distributions of spin objects which indi- 
cates the Fermi surface at different doping levels. In sec- 
tion we calculate the spin correlations which clearly 
demonstrate that antiferromagnetic order persists upon 
doping, and in section IVIIII we attempt to search for 
other exotic order when the antiferromagnetic correlation 
is weakened. Finally we give our conclusion in section llXl 



II. HOLE- AND ELECTRON-DOPED MODELS 

We start from the extended t-J model which was orig- 
inally proposed to describe hole-doped materials. 
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The nearest neighbor (nn) spin exchange constant J is 
fixed at 0.3. Farther than nearest neighbor hopping terms 
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are necessary in order to distinguish between hole and 
electron doping, tij = t,t', and t" when (ij) is a pair 
of sites at distances 1, v^, and 2 apart respectively, and 
is zero otherwise. The best fitting to ARPES results on 
SrsCuOaCla yields t = 1, t' = -0.3, and t" = 0.2ii 
In the case of hole doping, is a spin (or electron) 
creation operator . To describe electron-doped materi- 
als it is usual to apply the electron-hole transformation,"^ 
Cia — > aj^, where is a hole creation operator. The re- 
sulting Hamiltonian for electron-doped materials is iden- 
tical to but with Cia- replaced by a^g-, etc, and re- 
placed by —tij. As a result we can turn the Hamiltonian 
([T|) into an electron-doped model by flipping the signs of 
the hopping terms tij 3 Despite this similarity, one should 
be reminded that in the case of hole doping the operators 
Cia in the Hamiltonian are electron operators and the vac- 
uum state is a state with no electron. The condition of no 
double occupancy means that no more than one electron 
can occupy the same site. At half-filling, each site has 
exactly one electron and doping with holes means cre- 
ating vacancies by removing electrons. Translating into 
the language of the electron-doped model, the operators 
Cia are hole operators and the vacuum state is a state 
with no hole, i.e., it cannot accommodate any more elec- 
tron. The condition of no double occupancy means that 
each site can have no more than one hole. At half-filling, 
each site has exactly one hole and doping means filling 
up holes with electrons. To avoid confusion, we use the 
terms "spin objects" and "charge carriers" to describe ob- 
jects in the Hamiltonian In the case of hole doping, 
spin objects refer to electrons and charge carriers refer to 
holes. In the case of electron doping their meanings are 
reversed — spin objects refer to holes and charge carri- 
ers refer to electrons. In this paper, by "electron-doped 
model" we mean ^ with hopping parameters f = — 1, 
t' = 0.3, and r = -0.2. jl8j In principle hole-doped 
model should refer to the same extended t-J model with 
t = 1, t' = —0.3, and t" = 0.2. But due to complications 
caused by excited states of that modeljM we choose to 
compare the electron-doped model with the simple t-J 
model, i.e., with t = 1 and t' = t" = 0. We remark 
that the t-J model was originally proposed to describe 
hole-doped materials. But since they have different hop- 
ping terms, it will be unfair to conduct a quantitative 
comparison between the t-J and electron-doped models. 
Instead we will mostly concentrate on their qualitative 
differences. 



The electron-doped model with Hamiltonian Q and 
parameters J = 0.3, t — —1, t' = 0.3, and t" = —0.2 is 
solved by exact diagonalization on a square lattice with 
32 sites and periodic boundary conditions. TablelJshows 
the ground state energies and symmetries of the electron- 
doped model with one, two, and four electrons. Calcu- 
lations on the model with four electrons were performed 
on a cluster of AMD Opteron servers with 64 CPUs. 



TABLE I: Ground state energies, momenta and point group 
symmetries of the electron-doped model with Nc charge car- 
riers. Nb is the number of basis in that particular subspace. 
The ground state energy at half-filling is -11.329720. 





Nb 




k symmetry 


1 


150,297,603 


-13.913616 


(7r,0),(0,7r) 


2 


150,295,402 


-16.601689 


(0,0) d^2_y2 


4 


2,817,694,064 


-20.461647 


(0,0) s 



III. QUASI-PARTICLE DISPERSION IN THE 
ONE-ELECTRON SYSTEM 

In the electron-doped model the spectral function of 
spin objects at half-filling is defined as^^ 

^(k,c.) = - El + Et), (2) 

n 

where Eq and "i/jg are the ground state energy and wave 
function of the model at half-filling, and E^ and ij^^ are 
energy and wave function of the nth. excited state of the 
model with one charge carrier. v4(k, oj) can be evaluated 
easily using the continued fraction expansion. At each k 
point we use 300 iterations and add an artificial broad- 
ening factor e = 0.05 to the delta function. Fig. ^ shows 
A(k, lS) along three branches in the first Brillouin zone. 
At the Brillouin zone center (0,0), the spectral function 
has a broad structure. It does not have any noticeable 
low energy peaks. As we move along the (1,1) direc- 
tion towards (tt, tt) [Fig.^a)] the spectral weight spreads 
to lower and higher energies, forming well-defined peaks 
at low energies. When we go pass (7r/2,7r/2), the re- 
verse occurs and at (tt, tt), the spectral function has a 
broad structure again. A similar trend is observed in 
the branch from (0, 0) to (vr, tt) through (tt, 0) [Fig.IHb)]. 
This branch has the largest dispersion. A different trend 
is observed along the AFBZ boundary, (tt, 0) to (0, tt) 
[Fig.^c)]. There the spectral weight mostly concentrates 
in low energy states. As usual we define the quasiparticle 
weight as 

^ _ m\~C^a\r,)? .3. 
' «|£Uk.|V'o°)' 

where is the lowest energy one-electron state which 
has non-zero overlap with Cko-jV'o)- These values are tab- 
ulated in Table im together with the quasiparticle energy, 

£;(k) ^El~ El (4) 

which is the energy of the state -0^ in Eq.© relative to 
the ground state energy at half-filling. Most Zk are too 
small to make their corresponding peaks visible in Fig.^ 
Their positions are indicated by shaded arrows. An ob- 
vious exception is k = (tt, 0) which is the ground state 
momentum of the one-electron system. There the low- 
est energy peak has more than 60% of the total spectral 
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TABLE 11: Quasiparticle energy and weight in the electron- 
doped model. 



k 






(0,0) 


-1.452767 


0.00008 


/ 7r TT \ 
(4' 4,) 


-1.689698 


0.00354 


i TT TT \ 

\2' 2' 
(3tt Ztt\ 


-1.596649 


0.00465 


-1.669586 


0.00045 


(tt.tt) 


-1.450790 


0.00005 




-1.986268 


0.00237 


(^,0) 


-2.583895 


0.63608 


(f,0) 


-2.007321 


0.02755 


I 4 ' 4 / 


-2.092922 


0.01890 



weight. The quasiparticle energy is plotted in Fig. 12 The 
bandwidth is 3.777J, which is in qualitative agreement 
with the prediction of spin-polaron calculation^ 



IV. BINDING ENERGY 

The TVc-charge carrier binding energy is defined as the 
excess energy of a system with charge carriers over 
N(, single-carrier systems )Si 

Ei^^j ^ - £;0) - N,{El - (5) 

It indicates the tendency of the Nc charge carriers to 
form a bound state. From values in Tabled we find that 
E21 = -0.1042 and £"4/ = 1.2037. Compared to the t-J 
model where E21 = —0.0515,^^ it is tempting to interpret 
these numbers as evidence showing that two charge car- 
riers in the electron-doped model have larger tendency to 
form a bound pair than in the t-J model. However, we 
have two reasons to believe that this is not a fair conclu- 
sion. First of all we should not compare the tendency to 
form bound pairs in the two models based on the magni- 
tudes of their binding energies because they have different 
hopping terms. As we shall see in the next section, two 
charge carriers in the electron-doped model in fact have 
a smaller tendency to form a bound state than in the 
t-J model. Second, one must be careful in interpreting 
binding energies defined in Eq. JSJ because they are very 
susceptible to finite-size effects. Binding energies found 
in a finite system tend to be lower than their true values 
in the thermodynamic limit. As already pointed out in 
Ref. m we have no a priori reason to believe that E21 
can be extrapolated linearly in 1/iV, where N is the lat- 
tice size. Nevertheless, doing so with results at A'' = 16 
and 32 we obtain E21 ^ —0.01. This small value already 
hinted that the charge carriers may not form a bound 
state. When there are four charge carriers in the system, 
the large and positive En clearly shows that they have 
no tendency to form a bound state. 



V. CHARGE CARRIER CORRELATION IN 
REAL SPACE 

The real space correlation among charge carriers can be 
clearly displayed in the charge carrier correlation function 

C(r) = ((l-,v)(l-no)), (6) 

where the number operator of spin objects 

as in Eq. Note that we use the conventioni^ 

where Ne (r) is the number of equivalent pairs with sep- 
aration r. In this convention the correlation function 
satisfies the sum rule 

Y,NE{r)C{r)=N,~l, (8) 

and the probability for finding a pair of charge carriers 
at distance r apart is 

P(r) = NE{r)C(r)/{N, - 1). (9) 

Results in the two- and four-electron systems are shown 
m Fig. 13 Note that we use different symbols to dis- 
tinguish between two groups of correlations - those be- 
tween pairs of electrons in the same and opposite sub- 
lattices. A striking feature in the two-electron system 
is that the correlation between two electrons in oppo- 
site sublattices does not decay significantly with r. It is 
almost constant, implying that electrons in opposite sub- 
lattices are mostly uncorrelated. The correlation between 
two electrons in the same sublattice shows a very differ- 
ent trend. It is comparatively smaller than that in the 
other group and decays more significantly with distance 
r. The overall probability of finding a pair of electrons in 
the same and opposite sublattices are 0.3935 and 0.6065 
respectively. From these results we conclude that in the 
two-electron system, electrons prefer to stay in opposite 
sublattices where they can move almost independently 
of each other. Obviously this results from the fact that 
intra-sublattice hopping terms t' and t" in Ti. do not frus- 
trate the spin background and therefore allow electrons 
to move more freely. When we increase the number of 
electrons to four, the behaviors of the two groups of corre- 
lation become very similar. They show small fluctuations 
about the uncorrelated value {Nc — 1)/{N — 1), which is 
indicated by a dotted line in Fig. This shows that 
even in the four-electron system the electrons are mostly 
uncorrelated. The root-mean-square separation between 
two electrons \/{r^ are 2.2786 and 2.4131 in the two- 
and four-electron systems respectively. The proximity of 
these values to the root-mean-square separation between 
two uncorrelated electrons, 2.3827, again suggests that 
electrons in our systems are almost uncorrelated. 
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FIG. 1: (Color online) Spectral function at half-filling A(k,u} 
( 
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FIG. 2: Quasiparticle dispersion relation obtained from the 
spectral function at half-filling in Fig. The line is a guide 
to the eyes only. 



FIG. 3: (Color online) Charge carrier correlation in systems 
doped with two and four electrons. Squares and circles indi- 
cate that the pair of electrons are in the same and opposite 
sublattices respectively. Dotted lines are pair correlations of 
uncorrelated electrons in the respective cases. 



VI. CHARGE CARRIER CORRELATION IN 
MOMENTUM SPACE 

Next we go to the momentum space and study the 
momentum distribution function of spin objects (riko-) = 



(cj^^Cko-). One motivation for studying (nko-) is to learn 
about the Fermi surface of the model. For this purpose it 
is important to realize that in t-J-like models the momen- 
tum distribution function has a dome shape in the first 
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Brillouin zone. This feature results from minimizing the 
kinetic energji^SiS and is not related to the actual Fermi 
surface of the model. Nevertheless, those k points along 
the AFBZ boundary are not affected by this kinematic ef- 
fect. Therefore it is possible to extract information on the 
Fermi surface from these k points. A second motivation 
for studying the momentum distribution function is to see 
whether the single-carrier ground state is relevant to the 
physics of the multiply- doped system. This has been a 
subject of discussion in the t-J modeli^2*2iiS^ Note that in 
the electron-doped model some authors choose to report 
the momentum distribution function of charge carriers, 
instead of (nko-)- These distribution 



<>Ck<T4, 



functions have the following properties: 

M + {uka) = n^'^" = (N, + N,)/N, (10) 
5](nk.) = N^, (11) 



(12) 



where N„- is the number of spin objects with spin a. In 
this section we will start with the single-electron momen- 
tum distribution functions. We will show that they arc 
qualitatively similar to those in the t-J model. Thus their 
features are generic to J-like models. We will then dis- 
cuss momentum distribution functions of the two- and 
four-electron systems. We will show that they can be 
constructed from the single-carrier result and that the 
Fermi surfaces are consistent with small pockets at single- 
carrier ground state momenta. 



A. One-electron system 



(0,71) (71,71) (0,71) (7t,7t) 

0.2822 - - - 0.5297 - - - 0.5308 0.3754 - - - 0.4991 - - - 0.4997 

IN I I \ I 

I 0.5079 0.5309 I I 0.4864 0.4995 1 

IN I I N I 

, , \ 

0.4728 0.5145 0.5304 0.4653 0.4867 0.4981 

I Nil \ I 

\ , , \ , 

I 0.4762 0.5134 I I 0.4610 0.4824 1 

I '\ I I '\ I 

0.4615--- 0.4799--- 0.5288 0.4519--- 0.4590--- 0.1661 



(0,0) 



(71,0) (0,0) 



(a) 



(b) 



FIG. 4: Momentum distribution functions (a) (nkt) and (b) 
(rikj) in the ground state of the one-electron system with 
one spin-down object removed and momentum (vr, 0). Due to 
symmetry, only one quadrant of the Brillouin zone is shown. 



Fig.^b)], and an "antidips" in (jikt) is found at a k point 
which is displaced from the dip by the antiferromagnctic 
momentum (tTjTt) [(0,7r) in Fig.^fa)]. From Fig.^Jb) we 
find that the depth of the dip ('T-(,r,Tr)i)~ ('^(7r,o)i) is 0.334. 
This is very close to o)/2 which is 0.318, indicating its 
close tie with the Fermi surface. Furthermore, along 
the edge of the AFBZ are very small. Therefore our 
data is consistent with a small Fermi surface, or carrier 
pocket, at (7r,0). This is to be expected in a lightly 
doped antiferromagnet.^^ Note that all features described 
so far are qualitatively the same in the hole- and electron- 
doped models. They are generic features resulting from 
the kinematic effect and the antiferromagnctic order of 
the spin background. They do not reflect the different 
physics of the hole- and electron-doped models. 



Let us begin with the one-electron ground state with 
Sz — 1/2 and momentum (7r,0). (Note that Sz refers to 
the total spin of spin objects.) The momentum distribu- 
tion functions are shown in Fig. ^ We immediately no- 
tice two very prominent features that exist in both (rikt) 
and {riki): (i) there exist very sharp minima at (7r,0) 
or (0, tt); (ii) besides these sharp minima, they are of a 
dome shape with a maximum around (tt, tt) and slopes 
down towards a minimum at (0,0). 

As discussed above, the dome shape is a generic feature 
that also exists in the t-J modeli^^ The only difference 
is that the locations of the maximum and minimum of 
the "dome" are interchanged compared to those in the 
t-J model. This is obviously due to the opposite signs of 
t in the two models. This dome-shape feature therefore 
does not represent the shape of the true Fermi surface. 
Note that (nk|) and (rikx) shift above and below the half- 
filled value of 1/2 respectively due to the restriction from 
Eq. (nk.) < n'^^''. 

Similarly, sharp minima found in Fig. ^ are also found 
in the momentum distribution functions of the t-J model 
with one hole. Just like in the t-J model, a "dip" in 
(rikx) is found at the ground state momentum [(7r,0) in 



B. Systems with two and four electrons 

When there is an even number of spin objects, (jikt) = 
(jiki) and we drop the spin variable a from (rtko-). Fig- El 
shows (rik) in systems doped with two and four elec- 
trons. Again we can identity the same "dome-shape" 
structure found in the one-electron system. It is a 
generic feature of the model and does not reflect the 
physics of the charge carriers. Further evidence for this 
comes from the "height" of the dome, which is defined as 
An = (7T.(7r tt)) — (n(o,o))- In the two-electron system it is 
0.114, which is roughly the same as An-f -|- An^ = 0.117 
in the one-electron system. And in the four-electron sys- 
tem it is 0.225, roughly twice of that in the two-electron 
system. These agreements show that the dome-shaped 
structures at different doping levels are due to the same 
effect. 

Another feature common to the one-, two- and four- 
electron systems is that their (rik) have very prominent 
dips at k = (tt, 0) and (0,7r). These dips resemble elec- 
tron pockets. This is certainly not a generic feature of 
t- J-like models because it is not found in the t-J modelii^ 
Note that (tt, 0) and (0, tt) are along the AFBZ boundary 
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(0,71) (71,71) (0,71) (71,71) 

0.1739 0.5289 0.5306 0.1724 0.5574 0.5598 TABLE III: Charge-carrier distribution function (rTk) in sys- 

I 11^^ 1 terns doped with two and four electrons. Numbers in brackets 

I 4952 5304 I I 4287 5604 1 ^''-^ corresponding results from the additive approximation 

; ' \ ' ; ; ' \ ' ; Eq. lESJ. Note that in Eq. ^ we have to average the one- 

0.4380 o'".5009 0.5289 0.3517 o"'.5011 0.5574 carrier results over all degenerate ground states. 

I Nil \ I 

k two-electron system four-electron system 



I 0.4378 0.4952 i i 0.3623 0.4287 i 

I Nil N I 

I 

0.4163--- 0.4380--- 0.1739 0.3345--- 0.3517--- 0.1724 



(0,0) 0.1150 (0.1179) 0.2280 (0.2357) 

(f ,f ) 0.0935 (0.0941) 0.2002 (0.1882) 

(5,|) 0.0304 (0.0301) 0.0614 (0.0602) 

(^,%) 0.0008 (0.0009) 0.0021 (0.0018) 

(7r,7r) 0.0006 (0.0008) 0.0027 (0.0015) 

(7r,f ) 0.0024 (0.0021) 0.0051 (0.0042) 

, , (7r,0) 0.3574 (0.3550) 0.3901 (0.7099) 

FIG. 5: Momentum distribution function (rik) in systems j-tt q qq^2 (0 0928) 2108 (0 1855) 

with (a) two and (b) four electrons. (^^^ q q^^^ (0^0362) 0T338 (0^0724) 



(0,0) (71,0) (0,0) (71,0) 

(a) (b) 



where (rik) is not disguised by the generic dome-shape 
feature. Therefore these pocket-hke features should re- 
flect the physics of the systems. The fact that pocket- 
like features are found at these doping levels immediately 
suggests the relevance of the singly-doped state to the 
multiply-doped one. If electrons doped into the parent 
system behave like weakly interacting fermions, then it is 
reasonable to expect that multiple-electron systems can 
be approximated by filling up the single-electron band. 
This should lead to electron-pockets at single-electron 
ground state momenta. To make this argument quanti- 
tative, we consider the charge-carrier distribution func- 
tion (rTk). From Eq. pOfl . (nu) can be considered as 
the suppression of (nt) from its maximum value upon 
doping. If multiple-electron states can be built up from 
single-electron states, we expect the suppressions in (rik) 
due to individual electrons doped into the system to be 
additive, 



(nk)2 ^ ("-kT>i + («ki>i, 

(nk)4 ^ 2((nkT)i + (^ki>i,) (13) 



where (rTk)^^ is the distribution function of a system with 
Nc electrons. Table IIIII shows that in the two-electron 
system the additive approximation is satisfied at all avail- 
able k points. In the four-electron system it works satis- 
factorily at most k points. Obvious exceptions are (tt, 0) 
and (37r/4, 7r/4). Note that at this doping level Eq. fl^ 
cannot hold at (tt, 0) without violating Eq. H1U|) . As a re- 
sult, instead of filling states at (tt, 0) some electrons will 
fill the next available low energy states which, according 
to Fig. 13 are at (37r/4, 7r/4). Therefore our results show 
that the additive approximation works in the electron- 
doped model at doping levels up to at least 0.125. Note 
that this conclusion is consistent with that in section llVl 
if the electrons are uncorrelated, we expect that multiply- 
doped states can be build up from the singly-doped state. 



VII. SPIN ORDER 

Our previous discussions have been based on the sce- 
nario that antiferromagnetic spin order is preserved upon 
doping. This is a direct consequence of the fact that 
the t' and t" hopping terms do not frustrate the spin 
background. In this section we provide evidence for 
the existence of antiferromagnetic correlation. Antiferro- 
magnetic spin order can be measured directly using the 
spin correlation function (So • S^)- Results are shown in 
Fig. EJa). At half- filling the system is known to possess 
long-range antiferromagnetic order— and its spin corre- 
lation is shown as a reference. We note that in the two- 
electron system the spin correlation is not much weaker 
than that at half-filling, and more importantly it does 
not show significant decay beyond r = ^/2. This indi- 
cates that strong antiferromagnetic spin order exists in 
the system. The same qualitative trend is also observed 
in the four-electron system. Although the spin correla- 
tion is inevitably weakened due to higher doping level, 
it does not decay significantly beyond r = \pl. Another 
way to display the same data is through the static struc- 
ture factor, 

5(k) =^e'i'--(So-S,). (14) 

r 

Fig. Elb) shows the structure factors in systems doped 
with two and four electrons. As the doping level in- 
creases, the height of the antiferromagnetic peak at (tt, tt) 
is reduced. But it still remains prominent and there is no 
sign of enhancement at any other k point. Our results 
therefore indicate that antiferromagnetic order persists 
at least up to a; = 0.125. Note that this doping level is 
close to the point where the antiferromagnetic phase ends 
in the phase diagram of NdCeCuO^i which is a; = 0.13. 

VIII. CHARGE CURRENT CORRELATION 

The existence of a staggered pattern in the charge 
current correlation function has been established in the 
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FIG. 6: (Color online) (a) Spin correlation function and (b) 
static structure factor of systems with two and four electrons. 
Results at half- filling are given for reference purpose. In (a), 
empty and filled symbols represent positive and negative cor- 
relations respectively. 

t-J modeliS2iSS It has been interpreted as a direct ev- 
idence for the staggered-flux phase in the mean-field 
picturen^ In this picture the hopping motions of charge 
carriers frustrate the antifcrromagnetic spin background 
and lead to staggered chiral spin correlation. Binding 
of charge carriers is mediated through the attraction be- 
tween charge carriers with opposite vorticity. Numeri- 
cally it has been found that staggered current pattern 
exists in the two-hole t-J model with dx2-y2 symmetry 
only when the holes are loosely bound. It docs not ex- 
ist in states with other symmetry, nor when J/t is so 
large that the holes are tightly bound. The vorticity and 
charge correlations are found to be proportional to each 
other. However, in the electron-doped model a very ba- 
sic ingredient of the above picture is missing, namely, the 
hopping motions of charge carriers are mostly unfrustrat- 
ing. The result is that antifcrromagnetic order is robust 
and charge carriers do not have significant correlation. 
Consequently it is very unlikely that staggered current 
correlation can exist. However, as antifcrromagnetic spin 
order is weakened at higher doping level, more subtle cor- 
relations may emerge;— We therefore calculate the cur- 
rent correlation in the two- and four-carrier ground states 
and see if there exists a systematic trend as doping level 
increases. 

Fig. shows the spatial variation of the current corre- 
lation function (jkijmn), where 

jki = i{clci - cjck), (15) 

for systems doped with two and four electrons. Another 
way to display the same data is to define the vorticity 
V'(r) of a square plaquette by summing up the current 
around it in the counterclockwise direction. The vorticity 
correlation Cyy(r) = {V{r)V{0))/x is shown in Fig. |S1 
Our result for the two-electron system leaves little doubt 



that there is no staggered pattern in either {jkijmn) or 
Cvv {''')■ In the four-electron system there is again no 
clear indication of a staggered pattern. We therefore con- 
clude that intra-sublattice hopping terms in the electron- 
doped model do not favor the formation of a staggered 
pattern in the current correlation. This agrees with a 
recent mean-field study on the electron-doped model.— 
Furthermore, in the mean-field picture the loss of vor- 
ticity correlation implies that charge carrier binding is 
not favored, which agrees with our results in section [V] 
Finally we remark that in our four-electron system, it 
seems like at short distances the current correlation is 
stronger and exhibit a staggered pattern. This is most 
obvious when we compare Fig.|SJa) and (b). This seems 
to suggest that at doping level x = 0.125 the system 
may be starting to develop some other order due to the 
weakening of antifcrromagnetic correlation. However, the 
range within which we observe the "right" correlation is 
too short and the correlation is too weak for us to de- 
cide whether it has any significance. Therefore we are 
not able to make any definite statement concerning this 
matter. 



IX. CONCLUSION 

We have solved the electron-doped model with one, 
two, and four charge carriers on a 32-site square lat- 
tice. Our results cover doping levels up to a; = 0.125. 
In the electron-doped model, intra-sublattice hoppings 
of charge carriers do not frustrate the spin background. 
Most of our results presented above can be understood 
as consequences of this fact. Since hopping motions are 
mostly unfrustrating, charge carriers can propagate more 
freely. This is reflected in the large quasi-particle band- 
width in the singly-doped system. In systems doped with 
two and four electrons, the charge carrier correlation 
function shows that electrons are uncorrelated. Again 
this is due to the fact that electrons can hop more freely 
in the same sublattice. There is no evidence of charge 
carriers forming a bound state. Unfrustrating hopping 
motions also mean that antifcrromagnetic correlation in 
the spin background is better preserved upon doping. 
This is clearly shown in the spin correlation function and 
static structure factor. It also shows up in the Fermi 
surface of the system. In the singly-doped system quasi- 
particle weights at k points along the AFBZ boundary 
are small except at the single-carrier ground state mo- 
menta, i.e., (tTjO) and its equivalent points. This re- 
sembles a Fermi surface consisting of small pockets at 
single-carrier ground state momenta, which is expected 
in a lightly doped antiferromagnet.^® These small pock- 
ets persist in our systems doped with two and four elec- 
trons and are clearly visible in their momentum distri- 
bution functions of spin objects. Furthermore, momen- 
tum distribution functions of our multiply-doped systems 
can be well approximated by adding up the singly-doped 
momentum distribution functions. This re-assures that 
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FIG. 7: Current correlation {jkijmn) /x in the ground state of the (a) two-electron, and (b) four-electron systems. The reference 
bond mn is indicated by a shaded line. On other bonds, arrows point along the directions of positive correlation and line widths 
are proportional to (jkijmn) /x. For reference purpose, numerical values of {jkijmn) / x are shown next to some of the bonds. 
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FIG. 8: Vorticity correlation Cvvif) in the ground state 
of the (a) two-electron and (b) four-electron systems. The 
reference plaquette is indicated by a cross inside it. Those 
plaquettes touching the reference one do not have well-defined 
vorticity correlation because current operators on overlapping 
bonds do not commute. 



netic phase persists up to x = O.lSi ARPES experiment 
on the same material shows that before the Fermi sur- 
face becomes a large one centered at (tTjTt), small pock- 
ets will start to appear at (7r/2,7r/2) (and its equivalent 
points) as those at (tt, 0) evolve.^ However, the quasi- 
particle energy E{k) at k = (7r/2,7r/2) in the electron- 
doped model is quite high (see Fig.EJ. Assuming that the 
quasi-particle dispersion relation does not change much 
on further doping, electrons doped into the system will 
fill lower energy states at other k points first. Therefore it 
is not surprising that we do not see pockets developing at 
(7r/2, 7r/2). A recent theoretical calculation predicts that 
pockets will start to appear at (7r/2,7r/2) at a; = 0.144 
when the Fermi level crosses a different band from that 
at lower doping levels iS. Therefore it is possible that our 
doping level is not large enough to observe the change in 
the Fermi surface as revealed by ARPES experiment. 



charge carriers are uncorrelated. 

Our results show that antiferromagnetic order in 
electron-doped model persists at least up to doping level 
X = 0.125. We find no clear evidence of other orders 
existing in our systems. This is consistent with the 
phase diagram of Nd2-a;Cea;Cu04 whose antiferromag- 
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